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ABSTRACT 

We report on five binary pulsars discovered in the Parkes multibeam Galactic plane survey. All of the pulsars 
are old, with characteristic ages 1-Ilxl0 9 yr, and have relatively small inferred magnetic fields, 5-90xl0 8 G. 
The orbital periods range from 1.3 to 15 days. As a group these objects differ from the usual low-mass binary 
pulsars (LMBPs): their spin periods of 9-88 ms are relatively long; their companion masses, 0.2-1. 1M Q , are, 
in at least some cases, suggestive of CO or more massive white dwarfs; and some of the orbital eccentricities, 
10~ 5 <; e <^ 0.002, are unexpectedly large. We argue that these observed characteristics reflect binary evolution 
that is significantly different from that of LMBPs. We also note that intermediate-mass binary pulsars apparently 
have a smaller scale-height than LMBPs. 

Subject headings: Binaries: general — pulsars: individual (PSR J1232-6501, PSR J1435-6100, 
PSR J1454-5846, PSR J1810-2005, PSR J1904+0412) 



1. INTRODUCTION 

Most of the ~ 40 binary pulsars known in the disk of the 
Galaxy are millisecond pulsars with weak magnetic fields (B ~ 
10 8 G), spin periods 2 < P < 15 ms, and in nearly circular or- 
bits with companions of mass 0.15 <^,m2< 1 0.4 M Q , presumably 
He white dwarfs (WDs), some of which have been detected op- 
tically. These are the low-mass binary pulsars (LMBPs), and 
their formation mechanism is well understood. After a neu- 
tron star spins down to long periods and a low-mass companion 
evolves off the main sequence, a long phase of stable mass- 
transfer ensues, during which the syste m may be dete ctable as 
a low-mass X-ray binary (LMXB; see Verbunt 1993 for a re- 
view). Eventually the orbit is circularized (Phinney 1992), the 
pulsar spins up , its magnetic field is somehow quenched (e.g., 



Romani 1990| ), and a long-lived 'recycled' radio millisecond 



pulsar emerges. Despite some uncertainties, it appears that the 
birth -rates of LMXBs and LMBPs are comparable (L orimer 
2000), and this evolutionary model successfully accounts for 
many properties of LMBPs. However, it should be noted that 
20% of millisecond pulsars are isolated, and it is not clear how 
they have lost their presumed past companions. 

A small but growing group of binary pulsars consists of ob- 
jects with 15 < P < 200 ms, intermediate-mass companions 
(ni2 p> 0.5 Mq, likely CO or heavier WDs), and orbital eccen- 
tricities in some cases much larger than their LMBP counter- 
parts. These are the intermediate-mass binary pulsars (IMBPs), 
and it is not entirely clear how they fit into the evolutionary 
scheme outlined above. It has been suggested that such sys- 
tems undergo a period of unstable mass-transfer and common- 
envelope (CE) evolution ( van den Heuvel 1994 ). IMBPs may 
have more in common with the evolution of high-mass sys- 
tems that spend part of their lives as high-mass X-ray binaries 



(HMXBs) and are progenitors to eccentric-orbit double-neutron 
star binaries, with the difference that they were not sufficiently 
massive for a second supernova to have occurred. 

The vast majority of millisecond pulsars known in the Galac- 
tic disk is located within 2 kpc of the Sun. This is due to the loss 
of sensitivity of most surveys at larger distances, particularly 
along the Galactic plane. To probe the Galaxy-wide distribu- 
tion of LMBPs and to learn more about rare species of pulsars 
it is therefore desirable to search the distant Galactic plane with 
improved sensitivity. 

The Par kes multibeam survey (^yne et al. 2000|; M anchester 
etal. 2001) covers a region of the inner Galactic plane (\b\ <5°, 
-100° < I < 50°) with sensitivity far surpassing that of previous 
pulsar surveys. The main aim of the survey is to find young and 
distant pulsars, but it retains good sensitivity to fast-spinning 
pulsars. A radio frequency of 1374 MHz is used, reducing dele- 
terious propagation effects that affect the detectability of distant 
pulsars at low lati tudes. So far, the survey has discovered mor e 
than 500 pulsars ( |Camilo et al. 2000b| ; [Manchester et al. 2000| ). 
includi ng binary (Lyne et al. 200C ; Kaspi et al. 2000) and 
young ( Camilo et al. 2000a[ ) pulsars"^ 

In this Letter we report the discovery of five short-period pul- 
sars in binary systems. They contribute significantly to our un- 
derstanding of binary pulsar evolution and demographics. 

2. OBSERVATIONS AND RESULTS 



The survey uses the 13-beam receiver system at the 64-m 
Parkes telescope in NSW, Australia. Radio noise at a central 
frequency of 1374 MHz and spanning 288 MHz in bandwidth 
is filtered in a 96 x 3-MHz filter bank spectrometer in each of 
two linear polarizations, in observations lasting 35 min. Signals 
from complementary polarizations are added, and the 96 volt- 
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ages for each beam are sampled every 250 (is, digitized, and 
written to magnetic tape for off-line analysis. The data are then 
searched for periodic and dispersed s ignals using standard tech- 
niques (e.g., Manchester et al. 1996| ). 

Pulsars J1435-6100, J1810-2005, J1454-5846, 
J 1232-6501, and J 1904+04 12 were first detected in data col- 
lected on 1997 May 26, August 26, 1998 January 22, 24, and 
August 12, respectively. Following confirming observations, 
PSR J18 10-2005 has been monitored in a series of timing ob- 
servations with the 76-m Lovell telescope at Jodrell Bank, UK, 
while the remaining pulsars have been observed at Parkes. 

At Parkes we record data from the central beam in a man- 
ner otherwise identical to the survey observations, while track- 
ing each pulsar for about 15min on each observing day, 
with the exception that since MJD 51630 we have observed 
PSR J1435-6100 with a 512 x 0.5-MHz filter bank and a sam- 
pling interval of 125 /is at a central frequency of 1390MHz. 
Data are collected on a few days about every two months, co- 
inciding with epochs during which survey observations are in 
progress. Also, PSR J1904+0412 was observed on a monthly 
basis with the 305-m Arecibo telescope, from 1999 October 
through 2000 July, using the Penn State Pulsar Machine, a 
128 x 0.0625 MHz filter bank with 80/is sampling at a cen- 
tral frequency of 1400 MHz. The data, time-tagged with the 
start-time of the observations, are de-dispersed and folded at 
the predicted topocentric pulsar period, forming pulse pro- 
files; and pulse times-of-arrival (TOAs) are measured by cross- 
correlating these profiles with high signal-to-noise ratio (S/N) 
templates (Fig. [l]), created from the addition of many profiles. 
Similar procedures are used at Jodrell Bank, with the difference 
that the data are de-dispersed and folded on-line; also, 32 x 3- 
MHz filter banks were used until MJD 51400, and 64 x 1-MHz 
filter banks have been used since. 
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FIG. 1. — Integrated pulse profiles for five pulsars at a frequency of 
1374 MHz. The time resolution of each profile is indicated by a horizontal 
bar. The profiles for PSRs J1232-6501, J1454-5846, and J1904+0412 are the 
template profiles used to obtain TOAs. That for J1435-6100, at 1390 MHz, 
is the template for the high-resolution data obtained since MJD 51630, while 
that for Jl 8 10-2005 has time resolution a factor of six better than the template 
and most of the data used to obtain its timing solution. 

We then use the TEMPO timing software 9 to determine ce- 

9t 



lestial coordinates, spin, and orbital parameters for the pul- 
sars. This is done by first converting the measured TOAs to the 
barycenter using initial estimate s of pulsar para meters and the 
DE200 solar-system ephemeris ( iitandish 1982| ), and by mini- 
mizing timing residuals with respect to the model parameters. 
The parameters thus obtained are listed in Table |J and the cor- 
responding residuals are displayed in Figure || as a function of 
date. 
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FIG. 2. — Post-fit timing residuals as a function of date for the five bi- 
nary pulsars. All orbits have been well sampled, and residuals as a function of 
binary phase are featureless. 



The average flux densities listed in Table [T] were estimated 
by converting the observed S/N to a scale calibrated using sta- 
ble flux densities known for a group of high-dispersion measure 
(DM) pulsars. See Manchester et al. (2001) for further details 
of search and timing procedures. 



3. DISCUSSION 

3.1. Evolution of the new systems 

All five of the newly discoveredpulsars have low inferred 
magnetic fields (B < 10 10 G; Table pi) when compared with the 
vast majority of known pulsars (see Fig. [|), and all are in cir- 
cular binary systems. These characteristics indicate that all of 
the pulsars have interacted with their companions in the past, 
and have been recycled to some extent. However, their peri- 
ods and period derivatives (and hence B) are larger than those 
of most millisecond pulsars, as indicated in the P—P diagram 
of Figure |3|: the spin parameters of PSR J 1435-6 100 place it 
marginally within the group of LMBPs at the lower left of the 
diagram, while those for the remaining four pulsars place them 
squarely amidst the IMBPs and double-neutron star systems. 



See http://pulsar.princeton.edu/tempc 
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FIG. 3. — Observed period derivative vs. period for the subset of pul- 
sars in the Galactic disk with small period derivatives (more than 1000 known 
pulsars lie above P > 10~ 16 ). Dots denote isolated pulsars, circles indicate 
LMBPs, squares represent IMBPs, and diamonds depict high-eccentricity, 
double-neutron star binaries (see text). Large dots represent the pulsars pre- 
sented in this paper, labeled by their partial names. Two lines of constant in- 
ferred magnetic field strength and a line of characteristic age equal to 10 10 yr 

are indicated. Pulsars spur lip via mass arrrptinn must rpsirlp to the right of the 

spin-up limit (see Arzoumanian, Cordes, & Wasserman 1999 for a discussion). 



FIG. 4. — Orbital eccentricity vs. period of binary pulsars in the disk_pf the 
Galaxy with measured eccentricities e < 0.1. Symbols are as in Figure hi with 
triangles denoting upper limits. Three IMBPs mentioned in the text are identi- 
fied in small type. The dotted lines should contain 95% of the eccentricities of 
the LMBPs (circles), according to the model of Phinney & Kulkarni (1994). 

Tauris, van den Heuvel, & Savonije (2000) then extended 
this work to intermediate-mass (2-6 M Q ) donor stars. Remark- 
ably, they find that for a certain range of initial orbital peri- 
ods, such close binaries can survive periods of super-Eddington 
mass transfer on sub-thermal (few Myr) timescales without ex- 
periencing a CE phase. Depending on initial donor mass and 
orbital period, low-mass systems like the three we have discov- 
ered may result (see their Figs. 2 and 4). 

How shall we choose between these two alternative scenar- 
ios (low- versus intermediate-mass original companions)? De- 
spite their present low-mass companions, the newly discovered 
systems are unlikely to be standard LMBPs, as already noted, 
because of their relatively large P and B (PSR J 1904+04 12 
also has too large an eccentricity; Fig. |J). The intermediate- 
mass donor branch of evolution is therefore more suitable to 
explain the new systems: the intermediate-mass systems tend to 
have shorter and less stable periods of accretion, often at much 
higher rates, leading to a natural explanation for the larger P, B, 
and (in at least some cases) eccentricities. With this evolution- 
ary path, there is no need for a long-lived X-ray accretion phase. 
These systems might therefore not be descendants of standard 
LMXBs, and should be accounted for separately in birth-rate 
calculations. What the X-ray progenitors of such systems look 
like is of course an interesting and unresolved question. 

3.1.2. High-mass systems: common-envelope and a puzzle 

As is clear from Figure |J, the eccentricity of 
PSR J 1454-5846 is much higher than predicted by the con- 
vective fluctuation-dissipation theory of Phinney (1992). The 
pulsar therefore has P, B, and eccentricity larger than expected 
for LMBPs, and ;«2 ~ 1.1 Mq. We thus confidently classify 
it as an IMBP with a presumed O-Ne-Mg WD companion. 
It is likely to have undergone CE evolution and spiraled-in 
to its present Pb = 12.4 d from an initial period of several 
hundred days, with a companion of original mass 5-7 M Q 

These estimates of ni2 assume mi = 1 .35 Mq ( Thorsett & Chakrabarty 1999 ) and (' = 60° ; it is unlikely that more than one of the three systems is sufficiently 



Using the companion masses to attempt a classification of the 
new systems yields results which are mostly inconsistent with 
those derived from the spin parameters: PSR J1435-6100 has 
ma ~ 1.1 M , decidedly not compatible with a LMBP; of the 
remaining four systems only PSR J 1454-5846 (mi ~ 1.1 Mq) 
appears to be an IMBP, while the other three have 0.2 <,mx 

0.3M Q 10 — on this basis they should be classified as LMBPs, 
but their periods and magnetic fields are significantly larger 
than those of any LMBPs with remotely comparable binary pe- 
riods. 

One further piece of useful information is provided by the 
orbital eccentricities. Phinney (1992) derived a relationship be- 
tween eccentricity and binary period for LMBPs with Pb ;> 2 d 
that is remarkably consistent with observations. One key ingre- 
dient of the theory is that mass transfer to the neutron star via 
Roche-lobe overflow be stable over the giant phase of evolution 
of the com panion star. The relationshi p need therefore not hold 
for IMBPs ( |Phinney & Kulkarni 1994 , and f° r mree of me nve 
IMBPs wi th measured eccentricities identified so far (Camilo 
etal. 1996; |Tauris & Savonije 1999| ; [Edwards & Bailes 2000[ ) 
it does not (see Fig. Q). 

3.1.1. Low-mass systems: non-standard evolution? 

Tauris & Savonije (1999) considered the detailed non- 
conservative evolution of close binary systems with 1-2 M 
donor stars and accreting neutron stars, refining the well- 
known correlation bet ween Pt, and W2 for LMBPs (J oss, Rap- 
paport, & Lewis 1987). The three new low-mass systems 
(PSRs J 1232-6501, Jl 8 10-2005, and J 1904+04 12) follow this 
relation, considering the uncertainties in m^. 

io T 



face-on so as to have a > 0.4 Mq CO WD companion. 
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( iDewi & Tauris 2000|; Tauris, van den Heuvel, & Savonije 
2000). Edwards & Bailes (2000) recently reported the dis- 
covery of PSR Jl 157-51 12, a system broadly comparable to 
PSR J1454-5846, albeit with ft = 3.5 d and possibly a some- 
what larger companion mass. 

The pulsar J1435-6100 is likely to have a massive (m.2 ;> 
1 Ms) O-Ne-Mg WD companion, like PSR J1454-5846. It 
must have started with a very large orbital period so as not to 
coalesce during the CE/spiral-in phase, and ended with ft = 
1.35 d, much smaller than P b = 12.4d for PSR J1454-5846. A 
difficulty with understanding PSR J1435-6100 lies in its spin 
parameters: they are closer to those of LMBPs than IMBPs 
(see Fig. |h. In other words, despite a presumed short-lived 
(~ 10 4 yr) mass transfer phase in a CE (and hence very little 
accretion), the pulsar's magnetic field was somehow quenched 
to a very low value (5 x 10 8 G), while it was spun up to a fast 
initial rate (Pi <J 9 ms). Compare its parameters with those of 
the IMBP B0655+64: ft = 1.3 vs. 1.0 d; m 2 - 1.1 vs. 0.8 M ; 
both with similar eccentricities, and likely products of CE evo- 
lution. While the orbital parameters are thus fairly similar, 
the spin parameters are the most different within IMBPs: both 
B and the present-day period of PSR B0655+64 are 23 times 
larger than those of PSR J1435-6100. The re cently discovered 
PSR J1757-5322 ( fcdwards & Bailes 200dh has spin param- 
eters virtually identical to those of PSR J1435-6100 (Fig. g) 
and orbital parameters also similar to those of PSR B0655+64. 
The reason behind such contrasting sets of parameters between 
PSRs J1435-6100/J1757-5322 and B0655+64 is a puzzle. 

3.2. The scale-height of IMBPs and LMBPs 

The preceding discussion suggests that classifying pulsars 
by present-day companion mass alone may not be particu- 
larly useful. We therefore define IMBPs as objects that once 
had intermediate-mass donor stars. While this is a model- 
dependent definition, operationally it applies to pulsar systems 
with 10 < P < 200 ms and e < 0.01. Among these systems 
most, but not all (e.g., PSR J 1232-6501), have m 2 £ 0.4 M Q , 
and 0.5 £ ft £ 15 d. 

It is notable that seven of the 12 presently known IMBPs 
(squares in Fig. g) have been discovered in recent low- or 
inter mediate-latitude surveys (this Letter and E dwards & Bailes 
2000). This is despite the greater effective volume searched 
with at least comparable sensiti vity to pulsars with P > , 10 ms 



in some 'all-sky ' surveys (e.g., Camilo, Nice, & Taylor 1996 
Lyne et al. 1998] ). We now address this curiosity. 



The median distance perpendicular from the Galactic plane 
for 23 known LMBPs is |zL = 0.4 kpc, and three systems have 
\z\ > 1.8 kpc flCamilo 1999h . For the group of 12 IMBPs, 
\z\ m = 0.2 kpc and the largest distance is |zj = 0.5 kpc (C amilo 
1 999; pdwards & Bailes 2000| ; this Letter). Despite selection 
effects affecting these determinations for both populations, it 
appears that IMBPs have a smaller scale-height than LMBPs. 
The maximum perpendicular distance that a pulsar born near 
the plane can reach is approximately proportional to the square 
of its initial perpendicular velocity. A scale-height for IMBPs 
that may be a factor of 2-4 smaller than for LMBPs requires 
a velocity for IMBPs a factor of <; 2 smaller than for LMBPs. 
This is plausible, considering that a typical LMBP progenitor 
is a 1 + 1.3 Mq system while an IMBP may descend from a 
4 + 1 .3 M© system. In summary, the recent flurry of IMBP dis- 
coveries may be due simply to the fact that recent efforts are 
surveying with significant sensitivity where IMBPs tend to re- 
side — along the Galactic plane. Similar distributions apply 
to X-ray binaries: HMXB s have smaller average velocity an d 
scale-height than LMXBs ( |van Paradijs & McClintock 1995] ). 

The newly discovered IMBPs are distant objects (3 < d < 
10 kpc), and were detected because they are relatively luminous 
pulsars (2 <; L1400 <J 30mJykpc 2 ; Table [lb 11 . Therefore they 
need not contribute greatly to the overall population of binary 
pulsars in the Galaxy. However, in order to determine conclu- 
sively the scale-height of IMBPs, and their incidence among 
binary pulsars, it is necessary to perform careful modeling of 
the recent high-frequency surveys, and to measure proper mo- 
tions where possible. 
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Table 1 

Parameters for five binary pulsars 



PSRJ1232-6501 PSR J1435-6100 PSR J1454-5846 PSR J1810-2005 PSR J1904+0412 



R. A. (J2000) 


12 32 17.840(5) 


14 35 20.2765(4) 


14 54 10.908(2) 


18 10 58.988(2) 


19 04 31.382(4) 


Decl. (J2000) 


-65 01 03.33(4) 


-61 00 57.956(6) 


-58 46 34.74(3) 


-20 05 08.3(6) 


+04 12 05.9(1) 


Period, P (ms) 


88.2819082341(3) 


9.347972210248(6) 


45.24877299802(9) 


32.82224432571(9) 


71.0948973807(3) 


Period derivative, P 


8.1(2) X 10- 19 


2.45(4) x 10- 20 


8.16(7) x 10 -19 


1.51(7) x 10-' 9 


1.1(3) x 10- 19 


Epoch (MJD) 


51270.0 


51270.0 


51300.0 


51200.0 


51450.0 


Orbital period, P/, (days) 


1.86327241(8) 


1.354885217(2) 


12.4230655(2) 


15.0120197(9) 


14.934263(2) 


Projected semi-major axis, x (1-s) 


1.61402(6) 


6.184023(4) 


26.52890(4) 


11.97791(8) 


9.6348(1) 


Eccentricity, e 


0.00011(8) 


0.000010(2) 


0.001898(3) 


0.000025(13) 


0.00022(2) 


Time of ascending node, r asc (MJD)" .... 


51269.98417(2) 


51270.6084449(5) 


51303.833(4) 


51198.92979(2) 


51449.45(25) 


Longitude of periastron, u) (deg) 


129(45) 


10(6) 


310.1(1) 


159(30) 


350(6) 


Span of timing data (MJD) 


50940-51856 


50939-51856 


50981-51856 


50757-51817 


51089-51865 


Weighted rms timing residual (ps) 


200 


14 


100 


430 


240 


Dispersion measure, DM (cirT 3 pc) 


239.4(5) 


113.7(6) 


116.0(2) 


240.2(3) 


185.9(7) 


Flux density at 1400 MHz, S 140 o (mjy) . . 


0.3 


0.2 


0.2 


1.1 


0.3 


Derived parameters b 


Galactic longitude, / (deg) 


300.9 


315.2 


318.3 


10.5 


38.0 


Galactic latitude, b (deg) 


-2.2 


-0.6 


0.4 


-0.6 


-1.1 


Surface magnetic field strength, B (10 s G) 


90 


5 


60 


20 


30 


Characteristic age, t c (10 9 yr) 


2 


6 


0.9 


3 


11 


Mass function, f\ (Mq) 


0.0013 


0.1383 


0.1299 


0.0082 


0.0043 


Companion mass, mi (Mq ) 


>0.14 


>0.90 


>0.87 


> 0.28 


>0.22 


Distance, d (kpc) 


10 


3.3 


3.3 


4.0 


4.0 


Distance from Galactic plane, |z| (kpc) . . 
Radio luminosity, L1400 (mjykpe ) 


0.4 


0.03 


0.02 


0.04 


0.08 


30 


2 


2 


18 


5 



NOTE. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. Figures in parentheses are 
twice the nominal TEMPO uncertainties in the least-significant digits quoted, obtained after scaling TOA uncertainties to ensure \i = 1- 

"Due to the large covariance between ui and time of periastron (To) in standard TEMPO fits for pulsars with e<l, the solutions for pst?q J1232— 6501, 
J1435— 6100, and J 1 8 1 0— 2005 were obtained using the ELL1 model, where T dsc (uj = 0) and (ecosw,esina;) are fit instead ( [Lange et al. 2001 ). In these cases e 
and ui (as well as Tq) can be derived. For the other two pulsars we used the standard (BT) binary model that fits for e, w, and Tq — which is listed instead of T asc . 

b The following formulae are used to derive some parameters: B = 3.2 x 10 19 (PP) 1 / 2 G; t c = P/(2P); and /1 = j; 3 (2ir/ Pb) 2 Tq = (rri2 smVp /(m\ +m2) 2 , where 
Tq = GMq/c 3 =4.925/js, m\. m ; arp the pulsar anH rnmpaninn masses, respectively, and i is the orbital inclination angle. ni2 is obtained from the mass 



function, with mj = 1.35 Mq (Thorsett & Chakrabarty 1999|) and i < 90°. The distances are calculated from the DMs with the Taylor & Cordes (1993) free- 



electron distribution model; \z\ = ds'm \b\ ; and L1400 = Sl400^ 2 



